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In this paper, we study a two-dimensional model for the growth of a molecular clusters in superfluid 
helium at low temperature. In the model, particles of diameter a follow random ballistic moves of 
length 5 — a — 256a. Upon attachment on the cluster surface, particles allow one- step relaxation to 
the nearest two-fold coordinated site. Average coordination numbers of particles show that in the 
presence of relaxation the screening for incoming particles modifies the microscopic structure of the 
cluster. These results are in qualitative agreement with experimental aggregation of methane in He 
droplets, in which large abundance of fully coordinated sites is observed. The average coordination 
number increases with ^, showing that screening is relevant when relaxation is present. As the cluster 
size increases, the corresponding structure clusters changes from a compact to a fractal, typical of 
ballistic and diffusion limited models, respectively. A scaling ansatz describing the crossover between 
the two regimes is proposed. 
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I. INTRODUCTION 



Experiments with atoms and molecules embedded in He droplets of up to 10^ He atoms have recently became 
available.^ Single molecules were found to rotate freely in He droplets and are unique probes of superfluid helium 
on a microscopic scale!^ He droplets have extensively been used as hosts for the growth and spectroscopic study 
of small atomic and molecular clusters.^ It was obtained that experiments in He droplets often supersede other 
known techniques, such as free nozzle beam expansion, matrix isolation, and cryogenic cooling cells. In particular. 
He droplets offer straightforward control over the cluster size. One of the most spectacular observations was that 
of the self assembly of chain clusters of polar linear molecule and unusual ring clusters,^ which demonstrated the 
deep role of aggregation kinetics on the resulting cluster structure. These observations spawned expectations that 
He droplets could be used as a host in the formation of unique mesoscopic clusters consisting of several hundred 
and larger number of particles, which are interesting objects on their own and may be of practical use. Using of He 
droplet technique for assembly and investigation of large clusters is emerging. Formation of metal clusters having up 
to several thousand atoms in He droplets have been documented.^ ^ Recent infrared spectroscopic studies of ammonisP 
and methand^ clusters in He droplets consisting of up to 10^ molecules, both suggested formation of compact clusters. 
These observations are surprising taking into account the very low temperature of 0.38 K in the droplets, which we 
thought would favor formation of very porous clusters, as previously reported in the bulk liquid He.l^ 

Previous studies identified a few general models of cluster growth^ ^ such as ballistic aggregation (BA)P^^ and 
diffusion limited aggregation (DLA),^^ which give rise to densely branched and fractal clusters, respectively. However, 
to our knowledge those structures have never been observed in free atomic or molecular clusters, which are usually 
obtained at rather high temperature, and the structures - even the number of particles in the clusters - are defined by 
thermodynamics,^^ rather than aggregation kinetics. Helium droplets seem to be an ideal medium to study atomic 
and molecular aggregation, which allows the realization of some important paradigms such as adding one particle at 
a time and avoiding coalescence of clusters. Molecules in superfiuid helium are expected to move ballistically, until 
they recombine or are scattered on some elementary excitation in the superfiuid, which are extremely rare at the 
temperature of the droplets of just 0.38 K.^Thus, the mean free path of the molecules in He droplets should be very 
long. However, it remains unclear what are superfiuidity infiuences on the kinetics of the cluster formation inside the 
droplets. 

In order to obtain more insight into the mechanism of low temperature aggregation in He droplets, we calculated 
the structure of the two dimensional clusters using aggregation models modified in some important features. In 
standard models it is usually assumed that incoming particles immobilize immediately upon the impact. However, 
atoms and molecules may retain the ability to relax even at very low temperature due to quantum mechanical zero 
point delocalization and tunneling. The relaxation is guided towards some local minimum energy structures, which is 
defined by the bonding directionality, like in the case of covalent and hydrogen bonding, or by the largest coordination 
number such as in the case of Van der Waals bonding. In order to study the possible effect of superfluidity on the 
cluster structure, we allowed particles moving randomly with different mean free path before attachment to the cluster. 

The paper layout is as follows. In section |n| the model and the computation details are described. In section [llll 



the simulation results are presented and discussed. Finally, some conclusions are drawn in section IV 



II. THE MODEL 

Morphology of the cl usters obtained via aggregation of identical particles has been extensively studied experimentalljttH 
and theoretically) "'^^ l -'^^ -M^ DLA is probably the most studied aggregation model.l^ According to the standard DLA 
model,^ particles are released, one at a time, far from a cluster and follow unbiased Brownian motion until they 
touch the cluster where they instantly immobilize. DLA model generates structures with complex multiscale and 
multifractal properties. Several generalizations of the DLA model were proposed,— including models in which the 
particles take long random steps. 1^^^ For example, MeakiiP^ replaced the Brownian motion of the DLA model by 
Levy flights with the step length / distributed according to p{l) ~ The Meakin's model exhibits a continuous 
transition from BA- to DLA-like clusters as the exponent a is varied from tho inflnity. More recently, FerreirsP^ 
studied an extended on-lattice DLA model, in which the particles execute long flights with flxed length at random 
directions. This procedure reduces the screening^ for incoming particles, which is a central feature of the DLA-like 
models. DLA-like structures were observed for small flight length while a transition between compact and DLA-like 
structures was observed for long flights. 

The relaxation of particles upon attachment to the cluster was recently incorporated into the DLA calculations.^^ 
The calculations lead fractal clusters with dense "arms" characterized by a high coordination numbers of the particles 
as compared to those in the original DLA. The relaxation revealed to be important in aggregation of nanoemulsion 
droplets^^ as demonstrated by measured structure factor S{q) peaked at high values of the wavenumber which is 
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typical for a close packed structure. 

In the present work, we report the results of the two-dimensional calculations, which are based on an algorithm 
which has been described previously!!^ Simulations started with a single particle stuck to the origin. Additional 
particles with diameter a are released one at a time and execute ballistic steps of size S at random directions. The 
particles were launched from a circle of radius Ri which is taken to be Ri = Rmax + 5^, where Rmax is the largest 
distance from the origin among all particles in the cluster. When the walker touches any particle of the cluster its 
trajectory is terminated and the particle relaxes to the nearest position having two or more bonds. However, if the 
distance between the particle and the cluster becomes larger than Rk = 1005Rmax/ci^ the particle is discarded and a 
new one is released from the launching circle. Since the asymptotic behavior is typically reached in large clusters of 
about 10^ — 10^ particles, an efficient algorithm for off-lattice simulation is required. 

We adopted the standard DLA off-lattice simulation strategies.!^ First, we allow the particles outside a circle of 
radius Ri + bS taking long random steps of the length Text provided that the particle remains outside of the launching 
circle. We adopted rext = max[r — {Ri + 5^),(^], where r is the distance of the particle from the origin. Also, the 
movement in large empty areas in the inner region which surrounds the cluster (r < Rmax) are very computer time 
consuming, especially for large aggregates. Thus, we used a strategy that allows the particles inside the launching 
circle to take long steps of length Tint if they do not cross any part of the aggregate. Finally, the determination of 
the contacts between the wandering particles and the cluster in the off-lattice simulations is another time consuming 
task which was performed following algorithms described elsewhere.^ 



III. RESULTS AND DISCUSSION 

In contrast to the extensive studies concerning the mobility of electrons and He+ ions in liquid 
the motion of neutral atoms and molecules in liquid helium has not deserved much attention. It was obtained that 
their mobility is deter mined by scattering on roton and phonon elementary excitations in ^He,^^ and by quasi-particles 
at Fermi level in ^jj^jHHII ^j^^ mean free path of He+ ions in superfluid ^He at 0.5 K was estimated to be about 
one micron.*^ In turn, mean free paths of the order of about or less than the distance between helium atoms were 
obtained close to or above the superfluid transition temperature (2.17 K) in ^He^^ or in nonsuperfluid ^jj^jHHH] 
The scattering cross section for He+ was estimated to be about 5 x 10~^^ cm^, which corresponds to an effective 
radius of about 12 A.^^ The large size of the ions is due to the strong polarization interaction of the charge and He 
atoms, resulting in the formation of the "snowball" In comparison, interaction of neutral species with liquid helium 
is due to weak van der Waals forces, which primarily affects the distribution of He atoms in the first solvation shell.!^ 
Therefore the scattering cross section for atoms and small molecules, such as methane should be of the order or less 
than about 10~^^ cm^. Thus the mean free path of molecules in liquid ^He at T < 0.5 K should be about or greater 
than 10 /im, which is much larger than the particle diameter. Hence, particles experience long ballistic flight before 
collision with the next elementary excitation. In contrast, the motion of atoms and molecules at temperatures close 
to the superfluid transition, as well as in normal fluid ^He and ^He, should be diffusive with characteristic mean free 
paths of about or less than the particle's diameter. In our simulations, these experimental estimates correspond to 
the parameter 6 varying approximately from 1 up to 10^. 

Figure [l] shows the clusters grown for some values of the step size 6. The simulations were stopped when the 
cluster reached the border of a square region having size of 500a x 500a. As the parameter S increases, a continuous 
morphological transition from ramified to homogeneous clusters is observed. For large values of 5 the patterns are 
essentially those observed for the ballistic aggregation model but they become very similar to the DLA as 6 decreases. 
For comparison, clusters obtained via standard DLA and BA calculations without relaxation are shown in Fig. [ije) 
andjljf), respectively. It is seen that relaxation causes thicker branches of the aggregates which widths seem to be 
dependent on the S value. Inserts to panels in Fig. 1 show some examples of the local structure of the particles in the 
clusters. It is seen that relaxation favors large regions of close packing with some vacancies, contrasting with ordinary 
BA and DLA calculations which result in very erratic packing. 

Figure [2] shows the distribution of the number of the nearest neighbors, p^^ in clusters consisting of 10^ and 10^ 
particles in panels (a) and (b), respectively. For sake of comparison. Fig. [2|b) also includes the Pz distributions for 
the original DLA and BA calculations. It is seen that clusters obtained with relaxation have a large abundance of 
the fully coordinated {z = 6) and single vacancy sites {z = 5). In addition, sites having z = 3 and 4 correspond to 
the particles on the outer or inner boundary surfaces. The results clearly differ from the clusters obtained via BA 
and DLA algorithms both having very similar pz distributions peaking at z = 2. This indicates that the screening for 
incoming particles has no effect on the local arrangement of the particles in these models. Insert to Fig. |2|a) shows 
the probability of the occurrence of the fully coordinated site as a function of S. For each cluster size, the results 
indicate a logarithmic increase of the coordination number for small J, followed by saturation at large S. This is 
consistent with the BA-like scaling regime that is expected to be dominant for 5 larger than the characteristic cluster 
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size. These results are in qualitative agreement with experiments on the aggregation in He droplets, which show 
large abundance of fully coordinated sites. The overall shape of the distribution is shifted towards large coordination 
number as S increases, showing that screening becomes important when the relaxation mechanism is present. The 
insert to Fig. [2|a) also shows that the fraction of fully coordinated sites in the interior of the cluster remains nearly 
constant as more particles are added to the aggregate, in agreement with the experiments in He droplets showing 
compact clusters. 

The asymptotic DLA-like behavior is illustrated in figure 3, where sequential stages of the cluster grown with 
S = 128a are shown. At the early stages of the growth, the clusters closely resemble BA structures, which are 
characterized by homogeneous distribution of the branches. In larger clusters, the aggregates are dendritic resembling 
the DLA structures. A crossover from BA, at the early stages, to the DLA structures, at asymptotic stages, is observed 
at some critical cluster size which depends on the step size 5. 

An analysis of the cluster shape may provide some additional insights into the system. Figs. [I] and [s] suggest a 
morphological transition from nonfractal to fractal aggregates. Therefore, characteristics, such as fractality and/or 
multifractality,^^ as well as, the scaling exponents associated to the clusters^^ are important parameters describing 
overall cluster structure. 

In order to characterize the crossover from the BA- to the DLA-like structures, simulations were done on a square 
region having size of up to lO^a x lO^a with the step length ranging from S = la to 256a. The growth was stopped 
when the cluster reached the border of the limiting region. The fractal dimension, df, can be estimated using the 
mass-radius method by counting of the number of particles M enclosed by a circle of radius r centered at the origin 
taking the relation 



li df equals to the embedding dimension (of two in this work), the clusters are nonfractal. 

Figure |4ja) shows the mass-radius curves for several step sizes. For ^ > 1, the curves exhibit tenuous crossovers in 
the power laws, signaling the transition from the BA- to the DLA-like patterns. In all cases the value of df is close 
to 1.96(2) in small clusters and turns to df = 1.70(1) in larger clusters, which correspond to the fractal dimensions of 
the BA and the DLA aggregates, respectively.!^ Here and later in the paper, the uncertainties in the last significant 
digit are shown in parenthesis. The dependences on 6 of the mass M* and radius r* of the crossover cluster are shown 
in the insert to Fig. [4]^a). Scaling law relations, namely, r* ~ 5^ and M* ~ 5^ with z = 1.26(3) and a = 2.39(2), 
respectively, were found. 

A quantitative relation describing the crossover may be achieved using a scaling theory,^^ which establishes some a 
functional form for M(r, S) involving power laws. We propose a scaling function for the mass distribution: 



where f{x) ~ x^^, with df = dsA for x <C 1 and df = do la for x ^ 1. Here, dsA and do la are the fractal dimensions 
of the BA and DLA clusters. Eq. ([2| assumes that the crossover radius and mass scale as r* ~ 6^ and M* ~ 
respectively, and that the ratio r/r* rules the dependence of M with the cluster size. If the scaling hypothesis ([2| is 
correct, the plots Mj^^ against r/(5^ for all J values must collapse on a universal curve. 

The fractal dimensions used to determine the crossover points were dsA = 1-96 and d^LA = 1.7lj ^'^ l ^^ l The scaling 
hypothesis introduced by equation ([2| is corroborated in figure |4]^b), where all data collapsed on a universal curve are 
shown. Figure 4(a) shows that the mass of a cluster of size smaller than the crossover is independent of S implying the 
scaling relation a = zdsA- Indeed, the calculated values oi a = 2.39(2) and zdsA = 2.47(8) are in good agreement. 
The mass above of the crossover, r = 3000a, scales as M (5^ with = 0.26(3), obtained from the insert to Fig. 
[ij^b). Using the scaling behavior of the function ([2| and the exponents doLAi <^ and z from the fits, we obtained 
C = a — zd^LA = 0.24(8) which is in very good agreement with the value ^ = 0.26(3) above. Determining the scaling 
function and the corresponding power laws exponents, we have demonstrated the scale invariance of the clusters. 

The periphery or cluster border is defined as the particles forming an external layer impenetrable for incoming 
particles. Consequently, the spaces between adjacent particles at the border are smaller than a particle diameter. 
For the DLA clusters, almost all particles of the cluster belong to the border. Therefore, the number of peripheral 
particles, increases linearly with the cluster size. It is known that in the BA clusters S scales as 7V^-67(2) 34 rpj^^ 
calculated number of peripheral particles as a function of N for different 5 values is shown in figure [5] It is seen that 
upon increase of N , S exhibits a neat crossover from a scaling regime S ~ 7V^"^^ to S N . The exponent vba = 0.69 
is very close to that obtained in the original BA model, which corresponds to a growth faster than that observed in the 
compact and round Eden clusters having v = 0.5.1^ Collapses of all curves were obtained using a scaling hypothesis 
equivalent to Eq. ([2|: 



M(r) r' 



(1) 




(2) 
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where g{x) ~ with u = vba for x <C 1 and v = vdla = ^ ^oi x ^ 1. The scahng exponents used in the cohapse 
were a' = 1.62(4) and = 2.28(7). 



IV. CONCLUSIONS 

In this work, we modeled aggregation of spherical atomic or molecular particles in helium droplets. At temperatures 
of about 0.38 K the droplets are superfluid/ which presumably give rise to a long mean free path of the particles. We 
studied the dependence of the structure of the two dimensional clusters on the size of the random ballistic steps as 
well as on the relaxation of the particles in the cluster. The probability distribution of the local coordination number 
was used to measure the microscopic arrangement of the particles in the clusters. In the presence of relaxation the 
distribution was found to peak at the maximum coordination number of 2: = 6 contrasting with the standard BA and 
DLA models which peak at z = 2. Moreover, the fraction of particles with the largest coordination number increases 
with 6 at small and saturates for large values. 

Upon increase of the cluster size the morphology shows a crossover from nonfractal to fractal structures. The radial 
scaling and the size of the active growing zone were used to characterize this crossover. Scaling ansatz was proposed 
and successfully described the crossovers observed in the mass-radius and number of peripheral particles calculations, 
demonstrated by the collapses on universal curves observed for distinct values of the step length. The crossing point 
which delimits the fractal and non fractal regimes, increases with the mean free path. These results indicate that 
the structure of the clusters is a sensitive function of the mean free path of the molecules in superfluid He and that 
superfluid state of the clusters is essential for obtaining compact clusters as observed experimentally. 

In the future, we plan on expanding the calculations to 3D aggregation and modifying the standard DLA algorithm 
by taking into account the inter-particle interaction as well as exploring different relaxation mechanisms. Comparison 
of the experimental and theoretical results will give information on the kinetics of the motion of particles in liquid 
helium droplets, the role of inter-particle interactions as well as on relaxation pathways of the growing cluster at very 
low temperature. 
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FIGURE CAPTIONS 



FIG. 1. Clusters obtained with relaxation of the particles for different walker's step lengths: (a) ^ = a, (b) 6 = 4a, (c) 6 = 32a, 
and (d) 6 = 128a. Original DLA and BA models are shown in (e) and (f), respectively. The square region has a length 
L — 500a. The number of particles in the clusters is indicated in each panel. The insets show zooms of the regions indicated 
by the arrows. 



FIG. 2. (Color on-line) Distribution function of the coordination number in clusters consisting of N = 10^, and N — 10^ 
particles in panels (a) and (b), respectively, obtained at different step sizes 5. The results are averages over 100 simulation 
runs. The inset in (a) shows the fraction of particles having coordination number 2; = 6 as a function of 5 for two cluster sizes, 
N = 10^ (squares) and N = 10^ (circles). 



FIG. 3. A cluster at different stages of the growth for 5 — 128a. The number of particles in the clusters are shown in each 
panel. The boxes in (a) and (f) have lengths of 500a and lO^a, respectively. 



FIG. 4. (a) Mass versus radius for ^ = a, 4a, 16a, 64a and 256a from the bottom to the top. (b) Collapse of the data shown 
in panel (a). The inset in (a) shows the crossover radius r*, the mass for r = r*, and the power law fits (dashed lines). In (b), 
the inset shows a typical curve of the mass above the crossover point (r = 3000a) as a function of 5. The dashed line has a 
slope C = 0.30. Averages were done over 100 runs on a square region of size 5000a x 5000a in the main figures. 10 samples in 
a region of size 10^ a x 10^ a were used in order to obtain the insets. 



FIG. 5. Number of peripheral particles as a function of the total number of particles in the aggregate for 6 = a, 4a, 16a, 64a 
and 256a from the bottom to the top. The top inset shows the crossover total number of particles A/"* and the number of 
peripheral particles 5**. The bottom inset shows a typical curve for 5** above the crossover point (A^ = 2 x 10^) as a function 
of 5. Averages as in Fig. ^ 



